Critical local-moment fluctuations, anomalous exponents, and uj/T scaling in the 

Kondo problem with a pseudogap 
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Experiments in heavy-fermion metals and related theoretical work suggest that critical local-moment 
fluctuations can play an important role near a zero-temperature phase transition. We study such 
fluctuations at the quantum critical point of a Kondo impurity model in which the density of band 
states vanishes as |e| r at the Fermi energy (e = 0). The local spin response is described by a set 
of critical exponents that vary continuously with r. For < r < 1, the dynamical susceptibility 
exhibits uj/T scaling with a fractional exponent, implying that the critical point is interacting. 
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A number of stoichiometric (or nearly stoichiometric) 
heavy-fermion metals exhibit non-Fermi-liquid behavior 
when tuned to the vicinity of a magnetic quantum crit- 
ical point (QCP) jl|-f|. An important clue as to the 
nature of the quantum criticality has come from neu- 
tron scattering experiments [^|^] near the magnetic QCP 
of CeCue-zAuz. At the (rather small) critical Au con- 
centration, x c w 0.1, the dynamical spin susceptibil- 
ity is highly unusual in two respects: First, it satisfies 
lu/T scaling. Second, the frequency and temperature 
dependence obeys a fractional power law, described by 
the same anomalous exponent over essentially the en- 
tire Brillouin zone. There are indications that the stoi- 
chiometric system YbRh.2Si2 behaves similarly ||. These 
experiments directly suggest |^|| that the fluctuations 
of the individual local moments are also critical. While 
the standard Kondo behavior of local moments in simple 
metals has been studied extensively over the past four 
decades and is well understood Jl(j, the physics of crit- 
ical local-moment fluctuations is largely unexplored. It 
is therefore highly desirable to identify models that are 
amenable to controlled theoretical study. 

Further motivation for studying critical local-moment 
fluctuations comes from related theoretical work. We 
have recently shown pd| ] that competition between the 
Kondo and Ruderman-Kittel-Kasuya-Yosida interactions 
in a Kondo lattice model generates a new class of 
QCP, which we argue explains the aforementioned ex- 
periments in heavy fermions. Here, not only are the 
long-wavelength spin fluctuations critical but so also are 
the local-moment fluctuations; the weight of the Kondo 
resonance goes to zero at the QCP. The existence of criti- 
cal local fluctuations distinguishes such a "locally critical 
point" from the standard picture based on a spin-density- 
wave transition p^,[l3|]. To fully elucidate the proper- 
ties of the locally critical point, one must construct a 
Ginzburg-Landau description. This requires one to un- 
derstand the precise nature of the critical local mode 
that characterizes the destruction of the Kondo effect. 
The first step towards this goal is to develop an intu- 
ition about critical local-moment fluctuations in simpler 



models. 

This paper addresses just such a model: the single- 
impurity Kondo problem with a power-law pseudogap. 
The model has a quantum phase transition at a finite 
Kondo coupling We show that the QCP exhibits 

critical local-moment fluctuations and an associated de- 
struction of the Kondo effect very similar to those present 
at the locally critical point of the Kondo lattice; the local 
susceptibility displays co/T scaling with a fractional ex- 
ponent. The many-body spectrum of this model can be 
calculated exactly, an important virtue which should sig- 
nificantly aid the identification of the critical local mode. 

The Kondo model for a single spin-i impurity coupled 
to a conduction band is described by the Hamiltonian 
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where S is the impurity spin operator, cj^ creates an 
electron with spin z component a (= ±i) at the impurity 
site, and T l aa , (i — x, y, z) is the standard Pauli matrix. 
J is the Kondo exchange coupling, while V parametrizes 
nonmagnetic potential scattering from the impurity site. 

In the power-law version of this model, the conduction 
band is described by the (oversimplified) particle-hole- 
symmetric density of states pq] 
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for |e| < 1, 
for lei > 1. 



(2) 



In a metal (r = 0), any antiferromagnetic Kondo cou- 
pling J > causes the impurity moment to be quenched 
at temperature T = pffl. With a pseudogap (r > 0), 
by contrast, quenching occurs only for J > J c > Jli) . 
The strong-coupling (i.e., J > J c ) and weak-coupling 
(J < J c ) properties of the power-law Kondo model have 
been studied extensively p4] , p^| Here, instead, we 
study the critical behavior at J = J c (r,V) using nu- 
merical calculations and controlled analytical approxi- 
mations. (The model appears not to be integrable [Ml; 



1 



it also lacks conformal invariance.) This paper super- 
sedes an earlier preprint p2^| , which discussed only static 
critical properties, and focused on the large-iV limit of 
Eq. (EJ). More recent work ^3) has addressed QCPs in 
multi-channel Kondo problems with a pseudogap. 

Under conditions of strict particle-hole symmetry [V — 
in Eq. (EI)], a symmetric critical point (SCP) separates 
the weak- and strong-coupling regimes for all < r < i ; 
the strong-coupling regime and the SCP both vanish for 
r > h 10- The SCP is also encountered away from 
particle-hole symmetry for < r < r* sa 0.375; for r > 
r* , however, there is an asymmetric critical point (ACP) 
distinct from the SCP [Q. In all cases, the transition 
can be schematically represented as shown in Fig. EJ. 

Local vs impurity susceptibility. — Our analysis begins 
with the observation that the quantum critical behavior 
reveals itself, not in the response to a uniform magnetic 
field H, but rather in that to a local magnetic field h 
coupled solely to the impurity |24|,^5|]. These responses 
are measured, respectively, by the static impurity sus- 
ceptibility Ximp = — d 2 F imp /dH 2 \H=h=o, and the static 
local susceptibility xioc = —d 2 F- lrav / dh 2 \H=h=o 1 where 
-Pimp is the impurity contribution to the free energy. Nu- 
merical renormalization-group (NRG) results 0@ in- 
dicate that, whereas lim^^o TXim-p undergoes a jump as 
J passes through J c , liniT^o ^Xioc goes continuously to 
zero as the critical coupling is approached from below, 
and limr— »o ^Xioc = for all J > J c . (The same distinc- 
tion also holds in the large- iV limit p2[.) 

Static critical properties. — Given that the local field h 
(rather than the uniform field H) acts as a scaling vari- 
able, we define exponents (3, 7, <5, and x, describing the 
critical behavior of the the local susceptibility and the 
local-moment amplitude -Mi oc = {S z ) = — dF lmp /dh\H=o- 



M loc (J < J c , T = 0, h = 0) oc ( J c - Jf, 
Xloc(J> Jc,T=0) cx (J- J c )~ 7 , 
M loc (J = J c ,T = 0) cx H 1 / 5 , 

Xioc{J — Jc) — Cgtatic? 1 X - 
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Using a generalization JTq] of Wilson's NRG method 
p6| to treat the density of states in Eq. (|J), we have 
computed M\ oc and Xi oc for r between 0.1 and 2. For 
< r < 1, we find that M\ oc and xi oc obey Eqs. (||) near 
both the SCP and the ACP, establishing the continuous 
nature of these phase transitions 0]. For r > 1, by 
contrast, M\ oc undergoes a jump at the transition. 

Table Ej lists exponents for the SCP, along with their 
estimated nonsystematic (numerical-rounding and slope- 
fitting) errors. Data at J — J c exhibit power laws over at 
least five decades of h and T (e.g., see xi oc vs T in Fig. ||), 
allowing precise determination of 8 and x. The uncer- 
tainty in (3 and 7 is greater because rounding error cuts 
off the power laws as J approaches J c . Most runs were 
performed for an NRG discretization parameter A = 9, 
retaining all states within an energy 50T of the ground 
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FIG. 1. Schematic phase diagram showing the vicinity of 
the quantum critical point of the pseudogap Kondo model. 



state p5[ . To estimate the systematic discretization er- 
rors, a few runs were performed using a value A = 3 
lying closer to the continuum limit (A = 1) but requiring 
much more computer time. Critical exponents computed 
for A = 3 and A = 9 only narrowly fail to agree within 
their estimated nonsystematic errors, so we believe that 
the A = 9 exponents approximate the continuum values 
quite well. Restricting the power-law form of p(e) to a 
finite region around the Fermi energy to better approxi- 
mate real systems does not alter the critical exponents. 

The exponents listed in Table || have nontrivial r de- 
pendence. To better understand these exponents, we 
show that they satisfy certain hyperscaling relations, 
which can be derived in a standard fashion. We expect 
the singular component of the free energy to take the 
form |25fl 



^imp=T/(|J-J c |/T a , \h\/T» 



(4) 



Using Eq. (0), one readily finds that /? = (1 — b)/a, 7 = 
(26— 1) /a, 5 = b/ (1—b), and x — 26—1. These expressions 
lead to a pair of hyperscaling relations among the critical 
exponents, e.g., 

S=(l + x)/(l-x), (3 = j(l-x)/(2x). (5) 

In all cases, the exponents listed in Table EJ satisfy the 
hyperscaling relations to the accuracy of our calculations. 
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TABLE I. Properties of the symmetric critical point, ob- 
tained from NRG calculations. See the text for definitions 
of the exponents /3, 7, 8, and x (all determined using a 
band-discretization A = 9), and of the ratio R (calculated 
for A = 3). Parentheses surround the estimated nonsystem- 
atic error in the last digit (equal to 1 where omitted). 
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TABLE II. Exponents at the asymmetric critical point, 
from NRG calculations using a discretization parameter 
A = 9. The symbols are explained in Table @. 




i i i i i i u_ 

10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 10° 



T or cj 

FIG. 2. xioc(w = 0) vs T and Xk>c( T = 0) vs to at the 
symmetric critical point: NRG results for A = 3, r = 0.1 
and 0.4 Q. For a given r, xioc(^ = 0,T < 1) and 
Xk>c(M -C 1,T = 0) are described by equal exponents, x — y. 

Table || lists critical exponents at the ACP. Compar- 
ison with Table || shows that within their range of co- 
existence (0.375 < r < i), the SCP and ACP have dif- 
ferent exponents. For all r < 1, the hyperscaling rela- 
tions, Eqs. (||), are obeyed to within estimated errors. 
It proves difficult to determine the critical behavior for 
r = 1, where there are logarithmic corrections to scaling 
For 1 < r < 2, M\ oc is no longer critical, but xioc is 
described by exponents 7 = 2 — r and x = 1, values that 
are consistent with Eq. ([l]) if a = 1/(2 — r) and 6=1. 

Dynamical critical properties. — We have also com- 
puted the imaginary part of the dynamical local suscep- 
tibility, x'i oc (J = J c ,cj,T). Figure || shows some of our 
zero-temperature results at the SCP. The low-frequency 
NRG data at both the SCP and the ACP fit the form 

Xloc(J=Jc, W,T = 0) = CdynamicM^SgnW. (6) 

For 0<r< 1, we find y — x to within numerical error, an 
equality that is consistent with a scaling form: xioc(J — 
J C ,UJ,T) = T- x X(uj/T). Such an uj/T scaling cannot 
hold for 1 < r < 2, where we find y = 7 < x. 

For a given r between and 1, the data for Xk>c(^~ 
J c ,co,T) collapse onto a single function of to/T, as illus- 
trated in Fig. ||(a). This does not conclusively establish 
lo/T scaling because the NRG method is unreliable in 
the regime \lo\ ^ T. For small r, however, we can use a 
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FIG. 3. Comparison between NRG and small-r results for 
the symmetric critical point, (a) uo v X\ oc {J = J c) vs lo/T [ p5[ : 
NRG results for r = 0.4 and T/D = 10" 7 10" 6 , 10" 5 , and 
10 -4 (symbols) and the prediction of Eq. (g) (solid line), (b) 
Exponent 1 — x vs r: NRG data from Table M (symbols) and 
the prediction of Eq. ([)]), in the form of a spline fit (solid 
line) through numerical values for (poJ c ) 2 , extrapolated to 
the continuum limit A = 1 (see 

different approach to confirm the existence of scaling. 

For small r, the local spin-spin correlation function can 
be calculated algebraically by a procedure analogous to 
the standard e expansion [^7j : Since the critical coupling 
is small, the fixed point is accessible via an expansion 
in p J c m r. The unperturbed reference point (pqJ c = 
r = 0) describes the standard Kondo problem, so the 
perturbation series for xioc at the critical point contains 
logarithmic singularities. To leading logarithmic order, 
we find xioc(r) = |[1 — (poJc) 2 ln(7rTro/ sin7rTY)], where 
r « Pq. This leads to 

and to a dynamical spin susceptibility having the asymp- 
totic low-energy, low-frequency form 

„n Tn sin(7T?7/2) „ ft] to \ ,. 

XU^T)= i ( ^ T) Zj B(l- i —,l-r l ), (8) 

B being the Euler beta function. 

It follows from Eq. (||) that, for small r, the static local 
susceptibility and the imaginary part of the local suscep- 
tibility at T = have the forms specified by Eqs. (|^) 
and (0), respectively, with exponents 

x = y = 1-11 = I-(poJc) 2 , (9) 

and a universal (cutoff-independent) amplitude ratio 

C dynamic _ 7^-^(1-77/2) 
C static 2"sin( 7 r7 ? /2)r(l-7 7 )r(7 7 )r(7 7 /2)- 1 ' 
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These results can be compared with our NRG data: (i) 
As mentioned above, we find that y = x is obeyed to the 
accuracy of our calculations for all < r < 1. (ii) The 
calculated xioc obeys Eq. (SJ) within the range w » T 
where the NRG method is reliable [see Fig. (3 (a)], (iii) 
The numerical values of x agree remarkably well with the 
small-r result x = 1 — (poJ c ) 2 , even when r (and, hence, 
poJc) is not small [see Fig. 3(b)]. (iv) The R values in 
Table | fit Eq. (|l(]) to within 25%, a reasonable level of 
agreement given that the systematic errors in prefactors 
and critical couplings computed using the NRG are gen- 
erally greater than the errors in critical exponents. 

One of our key conclusions is that the dynamical spin 
susceptibility satisfies u)/T scaling, as shown by Eq. (||) 
for small r, and supported by the equality of the expo- 
nents x and y for all < r < 1. This has important 
implications for the field-theoretical description of the 
QCP. It indicates that a suitably defined relaxation rate 
is linear in temperature, which can only come about if the 
Ginzburg-Landau action, written in terms of the critical 
local modes, contains nonlinear couplings that are rele- 
vant in the renormalization-group sense |p8| . The fixed 
point must be interacting. For r > 1, the relaxation 
rate (oc T x / y ) is superlinear, which is consistent with a 
Gaussian fixed point. 

Finally, we note that the r — 1 pseudogap Kondo prob- 
lem is also important for impurities in high-T c supercon- 
ductors ^9j,^0|. The quantum critical regime will likely 
become accessible via scanning tunneling microscopy 
once measurements are extended to higher temperatures. 

In summary, we have combined numerical and analyt- 
ical approaches to obtain a consistent picture of the crit- 
ical properties of the Kondo problem with a conduction- 
electron density of states proportional to |e| r . At the 
QCP, the weight of the Kondo resonance has just gone 
to zero; as a result, the local-moment fluctuations are 
critical. In addition, the dynamical spin susceptibility at 
the QCP displays uj/T scaling with an anomalous expo- 
nent for all < r < 1. These features parallel those of the 
locally critical point in the Kondo lattice (TTJ . Thus, the 
pseudogap Kondo model provides a testing ground for 
studying critical local-moment fluctuations of the type 
seen in certain heavy- fermion metals The exact 

many-body spectrum of this impurity model can be cal- 
culated using NRG techniques, which should prove par- 
ticularly useful in identifying the proper local modes for 
characterizing the QCP. This, in turn, should shed much 
new light on the Ginzburg-Landau description of the lo- 
cally critical point of the Kondo lattice. These important 
issues are left for future work. 
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